\(\int \frac {x^4 (a+b \text {arccosh}(c x))^2}{(d-c^2 d x^2)^{5/2}} \, dx\) [215]

   Optimal result
   Rubi [A] (verified)
   Mathematica [A] (warning: unable to verify)
   Maple [B] (verified)
   Fricas [F]
   Sympy [F]
   Maxima [F]
   Giac [F]
   Mupad [F(-1)]

Optimal result

Integrand size = 29, antiderivative size = 482 \[ \int \frac {x^4 (a+b \text {arccosh}(c x))^2}{\left (d-c^2 d x^2\right )^{5/2}} \, dx=-\frac {b^2}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b^2 (1-c x)}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b^2 \sqrt {-1+c x} \sqrt {1+c x} \text {arccosh}(c x)}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b x^2 \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))}{3 c^3 d^2 \left (1-c^2 x^2\right ) \sqrt {d-c^2 d x^2}}+\frac {x^3 (a+b \text {arccosh}(c x))^2}{3 c^2 d \left (d-c^2 d x^2\right )^{3/2}}-\frac {x (a+b \text {arccosh}(c x))^2}{c^4 d^2 \sqrt {d-c^2 d x^2}}-\frac {4 \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))^2}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {\sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))^3}{3 b c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {8 b \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x)) \log \left (1-e^{2 \text {arccosh}(c x)}\right )}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {4 b^2 \sqrt {-1+c x} \sqrt {1+c x} \operatorname {PolyLog}\left (2,e^{2 \text {arccosh}(c x)}\right )}{3 c^5 d^2 \sqrt {d-c^2 d x^2}} \]

[Out]

1/3*x^3*(a+b*arccosh(c*x))^2/c^2/d/(-c^2*d*x^2+d)^(3/2)-1/3*b^2/c^5/d^2/(-c^2*d*x^2+d)^(1/2)+1/3*b^2*(-c*x+1)/
c^5/d^2/(-c^2*d*x^2+d)^(1/2)-x*(a+b*arccosh(c*x))^2/c^4/d^2/(-c^2*d*x^2+d)^(1/2)+1/3*b^2*arccosh(c*x)*(c*x-1)^
(1/2)*(c*x+1)^(1/2)/c^5/d^2/(-c^2*d*x^2+d)^(1/2)+1/3*b*x^2*(a+b*arccosh(c*x))*(c*x-1)^(1/2)*(c*x+1)^(1/2)/c^3/
d^2/(-c^2*x^2+1)/(-c^2*d*x^2+d)^(1/2)-4/3*(a+b*arccosh(c*x))^2*(c*x-1)^(1/2)*(c*x+1)^(1/2)/c^5/d^2/(-c^2*d*x^2
+d)^(1/2)+1/3*(a+b*arccosh(c*x))^3*(c*x-1)^(1/2)*(c*x+1)^(1/2)/b/c^5/d^2/(-c^2*d*x^2+d)^(1/2)+8/3*b*(a+b*arcco
sh(c*x))*ln(1-(c*x+(c*x-1)^(1/2)*(c*x+1)^(1/2))^2)*(c*x-1)^(1/2)*(c*x+1)^(1/2)/c^5/d^2/(-c^2*d*x^2+d)^(1/2)+4/
3*b^2*polylog(2,(c*x+(c*x-1)^(1/2)*(c*x+1)^(1/2))^2)*(c*x-1)^(1/2)*(c*x+1)^(1/2)/c^5/d^2/(-c^2*d*x^2+d)^(1/2)

Rubi [A] (verified)

Time = 0.67 (sec) , antiderivative size = 482, normalized size of antiderivative = 1.00, number of steps used = 20, number of rules used = 12, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.414, Rules used = {5934, 5892, 5912, 5913, 3797, 2221, 2317, 2438, 91, 12, 79, 54} \[ \int \frac {x^4 (a+b \text {arccosh}(c x))^2}{\left (d-c^2 d x^2\right )^{5/2}} \, dx=\frac {x^3 (a+b \text {arccosh}(c x))^2}{3 c^2 d \left (d-c^2 d x^2\right )^{3/2}}+\frac {\sqrt {c x-1} \sqrt {c x+1} (a+b \text {arccosh}(c x))^3}{3 b c^5 d^2 \sqrt {d-c^2 d x^2}}-\frac {4 \sqrt {c x-1} \sqrt {c x+1} (a+b \text {arccosh}(c x))^2}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {8 b \sqrt {c x-1} \sqrt {c x+1} \log \left (1-e^{2 \text {arccosh}(c x)}\right ) (a+b \text {arccosh}(c x))}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}-\frac {x (a+b \text {arccosh}(c x))^2}{c^4 d^2 \sqrt {d-c^2 d x^2}}+\frac {b x^2 \sqrt {c x-1} \sqrt {c x+1} (a+b \text {arccosh}(c x))}{3 c^3 d^2 \left (1-c^2 x^2\right ) \sqrt {d-c^2 d x^2}}+\frac {4 b^2 \sqrt {c x-1} \sqrt {c x+1} \operatorname {PolyLog}\left (2,e^{2 \text {arccosh}(c x)}\right )}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b^2 \sqrt {c x-1} \sqrt {c x+1} \text {arccosh}(c x)}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b^2 (1-c x)}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}-\frac {b^2}{3 c^5 d^2 \sqrt {d-c^2 d x^2}} \]

[In]

Int[(x^4*(a + b*ArcCosh[c*x])^2)/(d - c^2*d*x^2)^(5/2),x]

[Out]

-1/3*b^2/(c^5*d^2*Sqrt[d - c^2*d*x^2]) + (b^2*(1 - c*x))/(3*c^5*d^2*Sqrt[d - c^2*d*x^2]) + (b^2*Sqrt[-1 + c*x]
*Sqrt[1 + c*x]*ArcCosh[c*x])/(3*c^5*d^2*Sqrt[d - c^2*d*x^2]) + (b*x^2*Sqrt[-1 + c*x]*Sqrt[1 + c*x]*(a + b*ArcC
osh[c*x]))/(3*c^3*d^2*(1 - c^2*x^2)*Sqrt[d - c^2*d*x^2]) + (x^3*(a + b*ArcCosh[c*x])^2)/(3*c^2*d*(d - c^2*d*x^
2)^(3/2)) - (x*(a + b*ArcCosh[c*x])^2)/(c^4*d^2*Sqrt[d - c^2*d*x^2]) - (4*Sqrt[-1 + c*x]*Sqrt[1 + c*x]*(a + b*
ArcCosh[c*x])^2)/(3*c^5*d^2*Sqrt[d - c^2*d*x^2]) + (Sqrt[-1 + c*x]*Sqrt[1 + c*x]*(a + b*ArcCosh[c*x])^3)/(3*b*
c^5*d^2*Sqrt[d - c^2*d*x^2]) + (8*b*Sqrt[-1 + c*x]*Sqrt[1 + c*x]*(a + b*ArcCosh[c*x])*Log[1 - E^(2*ArcCosh[c*x
])])/(3*c^5*d^2*Sqrt[d - c^2*d*x^2]) + (4*b^2*Sqrt[-1 + c*x]*Sqrt[1 + c*x]*PolyLog[2, E^(2*ArcCosh[c*x])])/(3*
c^5*d^2*Sqrt[d - c^2*d*x^2])

Rule 12

Int[(a_)*(u_), x_Symbol] :> Dist[a, Int[u, x], x] /; FreeQ[a, x] &&  !MatchQ[u, (b_)*(v_) /; FreeQ[b, x]]

Rule 54

Int[1/(Sqrt[(a_) + (b_.)*(x_)]*Sqrt[(c_) + (d_.)*(x_)]), x_Symbol] :> Simp[ArcCosh[b*(x/a)]/b, x] /; FreeQ[{a,
 b, c, d}, x] && EqQ[a + c, 0] && EqQ[b - d, 0] && GtQ[a, 0]

Rule 79

Int[((a_.) + (b_.)*(x_))*((c_.) + (d_.)*(x_))^(n_.)*((e_.) + (f_.)*(x_))^(p_.), x_Symbol] :> Simp[(-(b*e - a*f
))*(c + d*x)^(n + 1)*((e + f*x)^(p + 1)/(f*(p + 1)*(c*f - d*e))), x] - Dist[(a*d*f*(n + p + 2) - b*(d*e*(n + 1
) + c*f*(p + 1)))/(f*(p + 1)*(c*f - d*e)), Int[(c + d*x)^n*(e + f*x)^(p + 1), x], x] /; FreeQ[{a, b, c, d, e,
f, n}, x] && LtQ[p, -1] && ( !LtQ[n, -1] || IntegerQ[p] ||  !(IntegerQ[n] ||  !(EqQ[e, 0] ||  !(EqQ[c, 0] || L
tQ[p, n]))))

Rule 91

Int[((a_.) + (b_.)*(x_))^2*((c_.) + (d_.)*(x_))^(n_.)*((e_.) + (f_.)*(x_))^(p_.), x_Symbol] :> Simp[(b*c - a*d
)^2*(c + d*x)^(n + 1)*((e + f*x)^(p + 1)/(d^2*(d*e - c*f)*(n + 1))), x] - Dist[1/(d^2*(d*e - c*f)*(n + 1)), In
t[(c + d*x)^(n + 1)*(e + f*x)^p*Simp[a^2*d^2*f*(n + p + 2) + b^2*c*(d*e*(n + 1) + c*f*(p + 1)) - 2*a*b*d*(d*e*
(n + 1) + c*f*(p + 1)) - b^2*d*(d*e - c*f)*(n + 1)*x, x], x], x] /; FreeQ[{a, b, c, d, e, f, n, p}, x] && (LtQ
[n, -1] || (EqQ[n + p + 3, 0] && NeQ[n, -1] && (SumSimplerQ[n, 1] ||  !SumSimplerQ[p, 1])))

Rule 2221

Int[(((F_)^((g_.)*((e_.) + (f_.)*(x_))))^(n_.)*((c_.) + (d_.)*(x_))^(m_.))/((a_) + (b_.)*((F_)^((g_.)*((e_.) +
 (f_.)*(x_))))^(n_.)), x_Symbol] :> Simp[((c + d*x)^m/(b*f*g*n*Log[F]))*Log[1 + b*((F^(g*(e + f*x)))^n/a)], x]
 - Dist[d*(m/(b*f*g*n*Log[F])), Int[(c + d*x)^(m - 1)*Log[1 + b*((F^(g*(e + f*x)))^n/a)], x], x] /; FreeQ[{F,
a, b, c, d, e, f, g, n}, x] && IGtQ[m, 0]

Rule 2317

Int[Log[(a_) + (b_.)*((F_)^((e_.)*((c_.) + (d_.)*(x_))))^(n_.)], x_Symbol] :> Dist[1/(d*e*n*Log[F]), Subst[Int
[Log[a + b*x]/x, x], x, (F^(e*(c + d*x)))^n], x] /; FreeQ[{F, a, b, c, d, e, n}, x] && GtQ[a, 0]

Rule 2438

Int[Log[(c_.)*((d_) + (e_.)*(x_)^(n_.))]/(x_), x_Symbol] :> Simp[-PolyLog[2, (-c)*e*x^n]/n, x] /; FreeQ[{c, d,
 e, n}, x] && EqQ[c*d, 1]

Rule 3797

Int[((c_.) + (d_.)*(x_))^(m_.)*tan[(e_.) + Pi*(k_.) + (Complex[0, fz_])*(f_.)*(x_)], x_Symbol] :> Simp[(-I)*((
c + d*x)^(m + 1)/(d*(m + 1))), x] + Dist[2*I, Int[((c + d*x)^m*(E^(2*((-I)*e + f*fz*x))/(1 + E^(2*((-I)*e + f*
fz*x))/E^(2*I*k*Pi))))/E^(2*I*k*Pi), x], x] /; FreeQ[{c, d, e, f, fz}, x] && IntegerQ[4*k] && IGtQ[m, 0]

Rule 5892

Int[((a_.) + ArcCosh[(c_.)*(x_)]*(b_.))^(n_.)/Sqrt[(d_) + (e_.)*(x_)^2], x_Symbol] :> Simp[(1/(b*c*(n + 1)))*S
imp[Sqrt[1 + c*x]*(Sqrt[-1 + c*x]/Sqrt[d + e*x^2])]*(a + b*ArcCosh[c*x])^(n + 1), x] /; FreeQ[{a, b, c, d, e,
n}, x] && EqQ[c^2*d + e, 0] && NeQ[n, -1]

Rule 5912

Int[((a_.) + ArcCosh[(c_.)*(x_)]*(b_.))^(n_.)*((f_.)*(x_))^(m_.)*((d1_) + (e1_.)*(x_))^(p_.)*((d2_) + (e2_.)*(
x_))^(p_.), x_Symbol] :> Int[(f*x)^m*(d1*d2 + e1*e2*x^2)^p*(a + b*ArcCosh[c*x])^n, x] /; FreeQ[{a, b, c, d1, e
1, d2, e2, f, m, n}, x] && EqQ[d2*e1 + d1*e2, 0] && IntegerQ[p]

Rule 5913

Int[(((a_.) + ArcCosh[(c_.)*(x_)]*(b_.))^(n_.)*(x_))/((d_) + (e_.)*(x_)^2), x_Symbol] :> Dist[1/e, Subst[Int[(
a + b*x)^n*Coth[x], x], x, ArcCosh[c*x]], x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[c^2*d + e, 0] && IGtQ[n, 0]

Rule 5934

Int[((a_.) + ArcCosh[(c_.)*(x_)]*(b_.))^(n_.)*((f_.)*(x_))^(m_)*((d_) + (e_.)*(x_)^2)^(p_), x_Symbol] :> Simp[
f*(f*x)^(m - 1)*(d + e*x^2)^(p + 1)*((a + b*ArcCosh[c*x])^n/(2*e*(p + 1))), x] + (-Dist[f^2*((m - 1)/(2*e*(p +
 1))), Int[(f*x)^(m - 2)*(d + e*x^2)^(p + 1)*(a + b*ArcCosh[c*x])^n, x], x] - Dist[b*f*(n/(2*c*(p + 1)))*Simp[
(d + e*x^2)^p/((1 + c*x)^p*(-1 + c*x)^p)], Int[(f*x)^(m - 1)*(1 + c*x)^(p + 1/2)*(-1 + c*x)^(p + 1/2)*(a + b*A
rcCosh[c*x])^(n - 1), x], x]) /; FreeQ[{a, b, c, d, e, f}, x] && EqQ[c^2*d + e, 0] && GtQ[n, 0] && LtQ[p, -1]
&& IGtQ[m, 1]

Rubi steps \begin{align*} \text {integral}& = \frac {x^3 (a+b \text {arccosh}(c x))^2}{3 c^2 d \left (d-c^2 d x^2\right )^{3/2}}-\frac {\int \frac {x^2 (a+b \text {arccosh}(c x))^2}{\left (d-c^2 d x^2\right )^{3/2}} \, dx}{c^2 d}+\frac {\left (2 b \sqrt {-1+c x} \sqrt {1+c x}\right ) \int \frac {x^3 (a+b \text {arccosh}(c x))}{(-1+c x)^2 (1+c x)^2} \, dx}{3 c d^2 \sqrt {d-c^2 d x^2}} \\ & = \frac {x^3 (a+b \text {arccosh}(c x))^2}{3 c^2 d \left (d-c^2 d x^2\right )^{3/2}}-\frac {x (a+b \text {arccosh}(c x))^2}{c^4 d^2 \sqrt {d-c^2 d x^2}}+\frac {\int \frac {(a+b \text {arccosh}(c x))^2}{\sqrt {d-c^2 d x^2}} \, dx}{c^4 d^2}+\frac {\left (2 b \sqrt {-1+c x} \sqrt {1+c x}\right ) \int \frac {x (a+b \text {arccosh}(c x))}{(-1+c x) (1+c x)} \, dx}{c^3 d^2 \sqrt {d-c^2 d x^2}}+\frac {\left (2 b \sqrt {-1+c x} \sqrt {1+c x}\right ) \int \frac {x^3 (a+b \text {arccosh}(c x))}{\left (-1+c^2 x^2\right )^2} \, dx}{3 c d^2 \sqrt {d-c^2 d x^2}} \\ & = \frac {b x^2 \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))}{3 c^3 d^2 \left (1-c^2 x^2\right ) \sqrt {d-c^2 d x^2}}+\frac {x^3 (a+b \text {arccosh}(c x))^2}{3 c^2 d \left (d-c^2 d x^2\right )^{3/2}}-\frac {x (a+b \text {arccosh}(c x))^2}{c^4 d^2 \sqrt {d-c^2 d x^2}}+\frac {\sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))^3}{3 b c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {\left (2 b \sqrt {-1+c x} \sqrt {1+c x}\right ) \int \frac {x (a+b \text {arccosh}(c x))}{-1+c^2 x^2} \, dx}{3 c^3 d^2 \sqrt {d-c^2 d x^2}}+\frac {\left (2 b \sqrt {-1+c x} \sqrt {1+c x}\right ) \int \frac {x (a+b \text {arccosh}(c x))}{-1+c^2 x^2} \, dx}{c^3 d^2 \sqrt {d-c^2 d x^2}}+\frac {\left (b^2 \sqrt {-1+c x} \sqrt {1+c x}\right ) \int \frac {x^2}{(-1+c x)^{3/2} (1+c x)^{3/2}} \, dx}{3 c^2 d^2 \sqrt {d-c^2 d x^2}} \\ & = -\frac {b^2}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b x^2 \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))}{3 c^3 d^2 \left (1-c^2 x^2\right ) \sqrt {d-c^2 d x^2}}+\frac {x^3 (a+b \text {arccosh}(c x))^2}{3 c^2 d \left (d-c^2 d x^2\right )^{3/2}}-\frac {x (a+b \text {arccosh}(c x))^2}{c^4 d^2 \sqrt {d-c^2 d x^2}}+\frac {\sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))^3}{3 b c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {\left (2 b \sqrt {-1+c x} \sqrt {1+c x}\right ) \text {Subst}(\int (a+b x) \coth (x) \, dx,x,\text {arccosh}(c x))}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {\left (2 b \sqrt {-1+c x} \sqrt {1+c x}\right ) \text {Subst}(\int (a+b x) \coth (x) \, dx,x,\text {arccosh}(c x))}{c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {\left (b^2 \sqrt {-1+c x} \sqrt {1+c x}\right ) \int \frac {c^2 x}{\sqrt {-1+c x} (1+c x)^{3/2}} \, dx}{3 c^5 d^2 \sqrt {d-c^2 d x^2}} \\ & = -\frac {b^2}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b x^2 \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))}{3 c^3 d^2 \left (1-c^2 x^2\right ) \sqrt {d-c^2 d x^2}}+\frac {x^3 (a+b \text {arccosh}(c x))^2}{3 c^2 d \left (d-c^2 d x^2\right )^{3/2}}-\frac {x (a+b \text {arccosh}(c x))^2}{c^4 d^2 \sqrt {d-c^2 d x^2}}-\frac {4 \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))^2}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {\sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))^3}{3 b c^5 d^2 \sqrt {d-c^2 d x^2}}-\frac {\left (4 b \sqrt {-1+c x} \sqrt {1+c x}\right ) \text {Subst}\left (\int \frac {e^{2 x} (a+b x)}{1-e^{2 x}} \, dx,x,\text {arccosh}(c x)\right )}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}-\frac {\left (4 b \sqrt {-1+c x} \sqrt {1+c x}\right ) \text {Subst}\left (\int \frac {e^{2 x} (a+b x)}{1-e^{2 x}} \, dx,x,\text {arccosh}(c x)\right )}{c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {\left (b^2 \sqrt {-1+c x} \sqrt {1+c x}\right ) \int \frac {x}{\sqrt {-1+c x} (1+c x)^{3/2}} \, dx}{3 c^3 d^2 \sqrt {d-c^2 d x^2}} \\ & = -\frac {b^2}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b^2 (1-c x)}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b x^2 \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))}{3 c^3 d^2 \left (1-c^2 x^2\right ) \sqrt {d-c^2 d x^2}}+\frac {x^3 (a+b \text {arccosh}(c x))^2}{3 c^2 d \left (d-c^2 d x^2\right )^{3/2}}-\frac {x (a+b \text {arccosh}(c x))^2}{c^4 d^2 \sqrt {d-c^2 d x^2}}-\frac {4 \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))^2}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {\sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))^3}{3 b c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {8 b \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x)) \log \left (1-e^{2 \text {arccosh}(c x)}\right )}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}-\frac {\left (2 b^2 \sqrt {-1+c x} \sqrt {1+c x}\right ) \text {Subst}\left (\int \log \left (1-e^{2 x}\right ) \, dx,x,\text {arccosh}(c x)\right )}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}-\frac {\left (2 b^2 \sqrt {-1+c x} \sqrt {1+c x}\right ) \text {Subst}\left (\int \log \left (1-e^{2 x}\right ) \, dx,x,\text {arccosh}(c x)\right )}{c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {\left (b^2 \sqrt {-1+c x} \sqrt {1+c x}\right ) \int \frac {1}{\sqrt {-1+c x} \sqrt {1+c x}} \, dx}{3 c^4 d^2 \sqrt {d-c^2 d x^2}} \\ & = -\frac {b^2}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b^2 (1-c x)}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b^2 \sqrt {-1+c x} \sqrt {1+c x} \text {arccosh}(c x)}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b x^2 \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))}{3 c^3 d^2 \left (1-c^2 x^2\right ) \sqrt {d-c^2 d x^2}}+\frac {x^3 (a+b \text {arccosh}(c x))^2}{3 c^2 d \left (d-c^2 d x^2\right )^{3/2}}-\frac {x (a+b \text {arccosh}(c x))^2}{c^4 d^2 \sqrt {d-c^2 d x^2}}-\frac {4 \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))^2}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {\sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))^3}{3 b c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {8 b \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x)) \log \left (1-e^{2 \text {arccosh}(c x)}\right )}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}-\frac {\left (b^2 \sqrt {-1+c x} \sqrt {1+c x}\right ) \text {Subst}\left (\int \frac {\log (1-x)}{x} \, dx,x,e^{2 \text {arccosh}(c x)}\right )}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}-\frac {\left (b^2 \sqrt {-1+c x} \sqrt {1+c x}\right ) \text {Subst}\left (\int \frac {\log (1-x)}{x} \, dx,x,e^{2 \text {arccosh}(c x)}\right )}{c^5 d^2 \sqrt {d-c^2 d x^2}} \\ & = -\frac {b^2}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b^2 (1-c x)}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b^2 \sqrt {-1+c x} \sqrt {1+c x} \text {arccosh}(c x)}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {b x^2 \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))}{3 c^3 d^2 \left (1-c^2 x^2\right ) \sqrt {d-c^2 d x^2}}+\frac {x^3 (a+b \text {arccosh}(c x))^2}{3 c^2 d \left (d-c^2 d x^2\right )^{3/2}}-\frac {x (a+b \text {arccosh}(c x))^2}{c^4 d^2 \sqrt {d-c^2 d x^2}}-\frac {4 \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))^2}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {\sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x))^3}{3 b c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {8 b \sqrt {-1+c x} \sqrt {1+c x} (a+b \text {arccosh}(c x)) \log \left (1-e^{2 \text {arccosh}(c x)}\right )}{3 c^5 d^2 \sqrt {d-c^2 d x^2}}+\frac {4 b^2 \sqrt {-1+c x} \sqrt {1+c x} \operatorname {PolyLog}\left (2,e^{2 \text {arccosh}(c x)}\right )}{3 c^5 d^2 \sqrt {d-c^2 d x^2}} \\ \end{align*}

Mathematica [A] (warning: unable to verify)

Time = 1.96 (sec) , antiderivative size = 382, normalized size of antiderivative = 0.79 \[ \int \frac {x^4 (a+b \text {arccosh}(c x))^2}{\left (d-c^2 d x^2\right )^{5/2}} \, dx=\frac {\frac {a^2 c x \left (-3+4 c^2 x^2\right ) \sqrt {d-c^2 d x^2}}{\left (-1+c^2 x^2\right )^2}-3 a^2 \sqrt {d} \arctan \left (\frac {c x \sqrt {d-c^2 d x^2}}{\sqrt {d} \left (-1+c^2 x^2\right )}\right )+\frac {a b d \left (-8 c x \text {arccosh}(c x)-\frac {\sqrt {\frac {-1+c x}{1+c x}} (1+c x)+2 c x \text {arccosh}(c x)}{-1+c^2 x^2}+\sqrt {\frac {-1+c x}{1+c x}} (1+c x) \left (3 \text {arccosh}(c x)^2+8 \log \left (\sqrt {\frac {-1+c x}{1+c x}} (1+c x)\right )\right )\right )}{\sqrt {d-c^2 d x^2}}+\frac {b^2 d \sqrt {\frac {-1+c x}{1+c x}} (1+c x) \left (-\frac {c x \left (-1+c^2 x^2+\left (-3+4 c^2 x^2\right ) \text {arccosh}(c x)^2\right )}{\left (\frac {-1+c x}{1+c x}\right )^{3/2} (1+c x)^3}+\text {arccosh}(c x) \left (\frac {1}{1-c^2 x^2}+\text {arccosh}(c x) (4+\text {arccosh}(c x))+8 \log \left (1-e^{-2 \text {arccosh}(c x)}\right )\right )-4 \operatorname {PolyLog}\left (2,e^{-2 \text {arccosh}(c x)}\right )\right )}{\sqrt {d-c^2 d x^2}}}{3 c^5 d^3} \]

[In]

Integrate[(x^4*(a + b*ArcCosh[c*x])^2)/(d - c^2*d*x^2)^(5/2),x]

[Out]

((a^2*c*x*(-3 + 4*c^2*x^2)*Sqrt[d - c^2*d*x^2])/(-1 + c^2*x^2)^2 - 3*a^2*Sqrt[d]*ArcTan[(c*x*Sqrt[d - c^2*d*x^
2])/(Sqrt[d]*(-1 + c^2*x^2))] + (a*b*d*(-8*c*x*ArcCosh[c*x] - (Sqrt[(-1 + c*x)/(1 + c*x)]*(1 + c*x) + 2*c*x*Ar
cCosh[c*x])/(-1 + c^2*x^2) + Sqrt[(-1 + c*x)/(1 + c*x)]*(1 + c*x)*(3*ArcCosh[c*x]^2 + 8*Log[Sqrt[(-1 + c*x)/(1
 + c*x)]*(1 + c*x)])))/Sqrt[d - c^2*d*x^2] + (b^2*d*Sqrt[(-1 + c*x)/(1 + c*x)]*(1 + c*x)*(-((c*x*(-1 + c^2*x^2
 + (-3 + 4*c^2*x^2)*ArcCosh[c*x]^2))/(((-1 + c*x)/(1 + c*x))^(3/2)*(1 + c*x)^3)) + ArcCosh[c*x]*((1 - c^2*x^2)
^(-1) + ArcCosh[c*x]*(4 + ArcCosh[c*x]) + 8*Log[1 - E^(-2*ArcCosh[c*x])]) - 4*PolyLog[2, E^(-2*ArcCosh[c*x])])
)/Sqrt[d - c^2*d*x^2])/(3*c^5*d^3)

Maple [B] (verified)

Leaf count of result is larger than twice the leaf count of optimal. \(2900\) vs. \(2(446)=892\).

Time = 1.35 (sec) , antiderivative size = 2901, normalized size of antiderivative = 6.02

method result size
default \(\text {Expression too large to display}\) \(2901\)
parts \(\text {Expression too large to display}\) \(2901\)

[In]

int(x^4*(a+b*arccosh(c*x))^2/(-c^2*d*x^2+d)^(5/2),x,method=_RETURNVERBOSE)

[Out]

4/3*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3*(c*x-1)*(c*x+1)*x^5+32*b^
2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3*c^2*arccosh(c*x)^2*x^7+16/3*b^2
*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3*c^2*arccosh(c*x)*x^7+181/3*b^2*(
-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3/c^2*arccosh(c*x)^2*x^3+40/3*b^2*(-
d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3/c^2*arccosh(c*x)*x^3-16*b^2*(-d*(c^
2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3/c^4*arccosh(c*x)^2*x-4*b^2*(-d*(c^2*x^2-
1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3/c^4*arccosh(c*x)*x+16/3*b^2*(-d*(c^2*x^2-1))^(
1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3/c^5*(c*x-1)^(1/2)*(c*x+1)^(1/2)-a^2/c^4/d^2*x/(-c^2
*d*x^2+d)^(1/2)+a^2/c^4/d^2/(c^2*d)^(1/2)*arctan((c^2*d)^(1/2)*x/(-c^2*d*x^2+d)^(1/2))+8/3*b^2*(-d*(c^2*x^2-1)
)^(1/2)*(c*x-1)^(1/2)*(c*x+1)^(1/2)/d^3/c^5/(c^2*x^2-1)*arccosh(c*x)^2-8/3*b^2*(-d*(c^2*x^2-1))^(1/2)*(c*x-1)^
(1/2)*(c*x+1)^(1/2)/d^3/c^5/(c^2*x^2-1)*polylog(2,c*x+(c*x-1)^(1/2)*(c*x+1)^(1/2))-8/3*b^2*(-d*(c^2*x^2-1))^(1
/2)*(c*x-1)^(1/2)*(c*x+1)^(1/2)/d^3/c^5/(c^2*x^2-1)*polylog(2,-c*x-(c*x-1)^(1/2)*(c*x+1)^(1/2))-16/3*b^2*(-d*(
c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3*(c*x-1)*(c*x+1)*arccosh(c*x)*x^5-55/3*
b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3/c^3*(c*x+1)^(1/2)*(c*x-1)^(1/
2)*x^2+64/3*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3/c^5*(c*x-1)^(1/2)
*(c*x+1)^(1/2)*arccosh(c*x)^2+4/3*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)
/d^3/c^4*(c*x-1)*(c*x+1)*x+16/3*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d
^3/c^5*(c*x-1)^(1/2)*(c*x+1)^(1/2)*arccosh(c*x)-8*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^
4-71*c^2*x^2+16)/d^3*c*(c*x-1)^(1/2)*(c*x+1)^(1/2)*x^6-8/3*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+1
18*c^4*x^4-71*c^2*x^2+16)/d^3/c^2*(c*x-1)*(c*x+1)*x^3+21*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118
*c^4*x^4-71*c^2*x^2+16)/d^3/c*(c*x-1)^(1/2)*(c*x+1)^(1/2)*x^4+1/3*a^2*x^3/c^2/d/(-c^2*d*x^2+d)^(3/2)-17*b^2*(-
d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3*x^5-4*b^2*(-d*(c^2*x^2-1))^(1/2)/(2
4*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3/c^4*(c*x-1)*(c*x+1)*arccosh(c*x)*x-13*b^2*(-d*(c^2*x^2-1))
^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3/c^3*(c*x-1)^(1/2)*(c*x+1)^(1/2)*arccosh(c*x)*x^2+
43/3*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3/c^2*x^3-4*b^2*(-d*(c^2*x
^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3/c^4*x-76*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8
*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3*arccosh(c*x)^2*x^5-44/3*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-
87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3*arccosh(c*x)*x^5+20/3*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*
x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3*c^2*x^7-1/3*a*b*(-d*(c^2*x^2-1))^(1/2)*(c*x-1)^(1/2)*(c*x+1)^(1/2)*(3*arcco
sh(c*x)^2*x^4*c^4-8*(c*x+1)^(1/2)*arccosh(c*x)*(c*x-1)^(1/2)*c^3*x^3-8*c^4*x^4*arccosh(c*x)+8*ln((c*x+(c*x-1)^
(1/2)*(c*x+1)^(1/2))^2-1)*x^4*c^4-6*arccosh(c*x)^2*x^2*c^2+6*(c*x+1)^(1/2)*arccosh(c*x)*(c*x-1)^(1/2)*c*x+16*c
^2*x^2*arccosh(c*x)-16*ln((c*x+(c*x-1)^(1/2)*(c*x+1)^(1/2))^2-1)*x^2*c^2-c^2*x^2+3*arccosh(c*x)^2-8*arccosh(c*
x)+8*ln((c*x+(c*x-1)^(1/2)*(c*x+1)^(1/2))^2-1)+1)/(c^6*x^6-3*c^4*x^4+3*c^2*x^2-1)/d^3/c^5-1/3*b^2*(-d*(c^2*x^2
-1))^(1/2)*(c*x-1)^(1/2)*(c*x+1)^(1/2)/d^3/c^5/(c^2*x^2-1)*arccosh(c*x)^3+84*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^
8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3/c*(c*x-1)^(1/2)*(c*x+1)^(1/2)*arccosh(c*x)^2*x^4+28/3*b^2*(-d*
(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3/c^2*(c*x-1)*(c*x+1)*arccosh(c*x)*x^3-
220/3*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3/c^3*(c*x-1)^(1/2)*(c*x+
1)^(1/2)*arccosh(c*x)^2*x^2+8*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x^4-71*c^2*x^2+16)/d^3
/c*(c*x-1)^(1/2)*(c*x+1)^(1/2)*arccosh(c*x)*x^4-32*b^2*(-d*(c^2*x^2-1))^(1/2)/(24*c^8*x^8-87*c^6*x^6+118*c^4*x
^4-71*c^2*x^2+16)/d^3*c*(c*x-1)^(1/2)*(c*x+1)^(1/2)*arccosh(c*x)^2*x^6-8/3*b^2*(-d*(c^2*x^2-1))^(1/2)*(c*x-1)^
(1/2)*(c*x+1)^(1/2)/d^3/c^5/(c^2*x^2-1)*arccosh(c*x)*ln(1+c*x+(c*x-1)^(1/2)*(c*x+1)^(1/2))-8/3*b^2*(-d*(c^2*x^
2-1))^(1/2)*(c*x-1)^(1/2)*(c*x+1)^(1/2)/d^3/c^5/(c^2*x^2-1)*arccosh(c*x)*ln(1-c*x-(c*x-1)^(1/2)*(c*x+1)^(1/2))

Fricas [F]

\[ \int \frac {x^4 (a+b \text {arccosh}(c x))^2}{\left (d-c^2 d x^2\right )^{5/2}} \, dx=\int { \frac {{\left (b \operatorname {arcosh}\left (c x\right ) + a\right )}^{2} x^{4}}{{\left (-c^{2} d x^{2} + d\right )}^{\frac {5}{2}}} \,d x } \]

[In]

integrate(x^4*(a+b*arccosh(c*x))^2/(-c^2*d*x^2+d)^(5/2),x, algorithm="fricas")

[Out]

integral(-(b^2*x^4*arccosh(c*x)^2 + 2*a*b*x^4*arccosh(c*x) + a^2*x^4)*sqrt(-c^2*d*x^2 + d)/(c^6*d^3*x^6 - 3*c^
4*d^3*x^4 + 3*c^2*d^3*x^2 - d^3), x)

Sympy [F]

\[ \int \frac {x^4 (a+b \text {arccosh}(c x))^2}{\left (d-c^2 d x^2\right )^{5/2}} \, dx=\int \frac {x^{4} \left (a + b \operatorname {acosh}{\left (c x \right )}\right )^{2}}{\left (- d \left (c x - 1\right ) \left (c x + 1\right )\right )^{\frac {5}{2}}}\, dx \]

[In]

integrate(x**4*(a+b*acosh(c*x))**2/(-c**2*d*x**2+d)**(5/2),x)

[Out]

Integral(x**4*(a + b*acosh(c*x))**2/(-d*(c*x - 1)*(c*x + 1))**(5/2), x)

Maxima [F]

\[ \int \frac {x^4 (a+b \text {arccosh}(c x))^2}{\left (d-c^2 d x^2\right )^{5/2}} \, dx=\int { \frac {{\left (b \operatorname {arcosh}\left (c x\right ) + a\right )}^{2} x^{4}}{{\left (-c^{2} d x^{2} + d\right )}^{\frac {5}{2}}} \,d x } \]

[In]

integrate(x^4*(a+b*arccosh(c*x))^2/(-c^2*d*x^2+d)^(5/2),x, algorithm="maxima")

[Out]

1/3*(x*(3*x^2/((-c^2*d*x^2 + d)^(3/2)*c^2*d) - 2/((-c^2*d*x^2 + d)^(3/2)*c^4*d)) - x/(sqrt(-c^2*d*x^2 + d)*c^4
*d^2) + 3*arcsin(c*x)/(c^5*d^(5/2)))*a^2 + integrate(b^2*x^4*log(c*x + sqrt(c*x + 1)*sqrt(c*x - 1))^2/(-c^2*d*
x^2 + d)^(5/2) + 2*a*b*x^4*log(c*x + sqrt(c*x + 1)*sqrt(c*x - 1))/(-c^2*d*x^2 + d)^(5/2), x)

Giac [F]

\[ \int \frac {x^4 (a+b \text {arccosh}(c x))^2}{\left (d-c^2 d x^2\right )^{5/2}} \, dx=\int { \frac {{\left (b \operatorname {arcosh}\left (c x\right ) + a\right )}^{2} x^{4}}{{\left (-c^{2} d x^{2} + d\right )}^{\frac {5}{2}}} \,d x } \]

[In]

integrate(x^4*(a+b*arccosh(c*x))^2/(-c^2*d*x^2+d)^(5/2),x, algorithm="giac")

[Out]

integrate((b*arccosh(c*x) + a)^2*x^4/(-c^2*d*x^2 + d)^(5/2), x)

Mupad [F(-1)]

Timed out. \[ \int \frac {x^4 (a+b \text {arccosh}(c x))^2}{\left (d-c^2 d x^2\right )^{5/2}} \, dx=\int \frac {x^4\,{\left (a+b\,\mathrm {acosh}\left (c\,x\right )\right )}^2}{{\left (d-c^2\,d\,x^2\right )}^{5/2}} \,d x \]

[In]

int((x^4*(a + b*acosh(c*x))^2)/(d - c^2*d*x^2)^(5/2),x)

[Out]

int((x^4*(a + b*acosh(c*x))^2)/(d - c^2*d*x^2)^(5/2), x)